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Abstract. The spin structures of small spin-1 condensates (N < 1000) under a 
magnetic field B has been studied beyond the mean field theory (MFT). Instead of 
the spinors, the many body spin-eigenstates have been obtained. We have defined 
and calculated the spin correlative probabilities to extract information from these 
eigenstates. The correlation coefficients and the fidelity susceptibility have also been 
calculated. Thereby the details of the spin-structures responding to the variation of 
B can be better understood. In particular, from the correlation coefficients which is 
the ratio of the 2-body probability to the product of two 1-body probabilities, strong 
correlation domains (SCD) of B are found. The emphasis is placed on the sensitivity of 
the condensates against B. No phase transitions in spin-structures are found. However, 
abrupt changes in the derivatives of observables (correlative probabilities ) are found 
in some particular domains of B. In these domains the condensates are highly sensitive 
to B. The effect of temperature is considered. The probabilities defined in the paper 
can work as a bridge to relate theories and experiments. Therefore, they can be used 
to discriminate various spin-structures and refine the interactions. 
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1. Introduction 

Bose-Einstein condensates are ideal artificial systems for quantum manipulation. The 
spinor condensates have very rich spin structures, and their swift response to the external 
field is notable pQ. Usually, a condensate would contain more than 10 4 atoms. Due to 
the progress in techniques, much smaller condensates (say, particle number N < 1000) 
could be produced. Note that the strength of the effective interaction between the 
particles C2 is proportional to the average particle density, which is proportional to 
jy-3/5 (evaluated based on the Thomas- Fermi approximation). Thus the particles are 
subjected to a stronger interaction in smaller systems. In fact, in the study of spin- 
dynamics, the rate of evolution depends on r = C^tjh, rather than t [2], where t is 
the time. Therefore, the evolution will be swifter in the smaller systems. Thus they 
might be even more suitable for manipulation. In principle, these systems could be 
more precisely prepared (such as iV and M, the total magnetization, could be more 
rigorously given). Therefore, some properties not contained in large condensates might 
emerge (say, if N < 50, whether N is even or odd might be serious), and some properties 
might depend on the external field very sensitively in some particular cases (as shown 
below) . 

There are already a number of literatures dedicated to the large condensates of 
spin-1 atoms. The commonly used theoretical tool is the mean field theory (MFT) 
[31 HI El EJ [71 [HI [9] . The ground states with the ferromagnetic phase and polar phase are 
found. Their modes of excitation and their dynamical behavior have been studied, and 
very rich physical phenomena have been found (say, the formation of spin-domains and 
spin- vortices). On the other hand, for small condensates, a theory goes beyond the MFT 
might be more appropriate. In this paper, the spin-structures of spin-1 small condensates 
under a magnetic field B are studied by using a many-body theory in which the parentage 
coefficients of spin-states are used as a tool [101 EE]. The variation of B is considered as 
adiabatic and the discussion is limited to static behavior. One-body and two-body spin 
correlative probabilities are defined and calculated. They are used to extract information 
from the spin eigenstates. [12] Furthermore, the correlation coefficients are defined and 
calculated to measure quantitatively the spin correlation, and the fidelity susceptibilities 
are also calculated to measure quantitatively the sensitivity of the ground states against 
the change of B. Thereby a detailed and deeper description on spin-spin correlation 
has been obtained that might lead to a better understanding on spin structures. In 
particular, the spin probabilities defined in this paper are observables, they might 
serve as a bridge to relate theories and experiments, therefore can be used to clarify 
various spin-structures and interactions. The emphasis is placed on the response of the 
condensates to B. The related knowledge might be useful for quantum manipulation. 
Both the cases with the temperature T zero and nonzero are considered. 
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2. Hamiltonian and the eigenstates 

Let N spin-1 atoms be confined by an isotropic and parabolic trap with frequency u. 
The interaction is Vy = 5{vi — rj)(c + c 2 f i ■ fj), where and f i are, respectively, 
the position vector and the spin operator of the z-th particle. A magnetic field B lying 
along the Z-axis is applied. We consider the case that the size of the condensate is 
small so that it is smaller than the spin healing length. In this case, the single mode 
approximation (SMA) is applicable. p3] Under the SMA all particles will have the same 
spatial wave function <p(r). After the integration over the spatial degrees of freedom, 
we arrive at a model Hamiltonian 

H mod = C 2 S 2 -pJ2fi* + Q E(/-) 2 > (!) 

i i 

where S is the total spin operator of the N particles, C 2 = c 2 J dr|0(r)| 4 /2, f iz is the 
^-component of f i: p = —^bB/2, q = (fisB) 2 /4Em, the Bohr magneton and Em 
the hyperfine splitting energy. The last two terms of H mod are the linear and quadratic 
Zeeman energies, respectively. 

The eigenstates of H mod will have 5* and its z-component M (namely, the 
magnetization, M > is assumed) conserved when the quadratic term is neglected. 
Therefore they can be denoted as $smi S = N, N — 2, to 1 or 0. It has been proved 
that ti$ M is unique without further degeneracy |14j . They together form a complete set 
for all the totally symmetric spin-states of / = 1 systems. When the quadratic term is 
taken into account, M will remain to be conserved, but S will not. In this case "&g M can 
be used as basis functions for the diagonalization of i? mo d- In this approach, a powerful 
tool, the fractional parentage coefficients, that we have developed previously is used for 
the calculation of related matrix elements [TOl [TT] . 

Using these coefficients, a particle (say, the particle 1) can be extracted from -&g M 

as 

where X/X-Q ^ s the spin-state of the particle 1 in component \i — or ±1. The fractional 
parentage coefficients have analytical forms as 



AN] 

a SMn 



[l + (-l)»-*}(N-S)(S + l) SM 
\ 2N(2S + 1) s 



Cfu M 5+l,M-u5 (3) 



b {N} 



[l + (-l)N-s] S (N + S + l) „ SM ( . 
\ 2N(2S+1) u i»,s-i,m-v W 

where the Clebsch-Gordan coefficients are introduced. 

Once a particle has been extracted, the calculation of the matrix elements is straight 
forward, and we have 

(^M'l#modfe> = hr,Ai{^s[C 2 S(S + 1) - P M] + qQ™}, (5) 
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where 



Qs's = N(^, M \(f lz f\^ M ) = ArSW- (6) 



NMfj. _ r r, {N} , 2 , , h {N} v 2 i , r {JV} i{N} 

QS'S - d S'S[{asM^) + \PSMn) \ + °S', S-2%- 2 , Mfj,°SMfi 

r {N} ,{AT} 1 /_s 

For an arbitrary iV, after a procedure of diagonalization, the set of eigenenergies and 
the corresponding eigenstates 0«Af = J2 S ds^^SM can obtained, where Ei is in the 
order of increasing energy. Since the set "dg M is complete, the set G^m is exact for H mod . 

3. Spin correlative probabilities 

For the condensates with nonzero spins, it has been shown theoretically that there 
are various spin-structures. In order to confirm these structures experimentally, one 
has to define some measurable physical quantities. For the spatial degrees of freedom, 
it is reminded that the one-body density p(r) can provide information on the spatial 
distribution of the particles, and the two-body density p(r±, r 2 ) can describe the spatial 
correlation. Similar quantities can be defined in the spin space. Each spin eigenstate 
can be written as 

e jM = E^M = Ex,(iH M (8) 

s fi 

where (p^ can be obtained via Equation (T5]). The normality 1 = (GWlGW) = 
^2 l j,( ( P t /^\ ( P t ^ d ) implies that (y^ M |<^ M ) is the probability of particle 1 in /i-component 
. Then, we define the one-body probability 

^ = W> = EW? ( 9 ) 

SS' 



Note that NP^" 1 is just the average population of the /i-component, and is an observable 
which can be directly measured via the Stern-Gerlach technique. 

In the case with B = 0, S is a good quantum number, and the i-th state has 
Si = N + 2 - 2z (if c 2 < 0), or # = 2{i - 1) (if c 2 > 0).[5J In this case, df 1 = Ssftdgf , 
and we have 

pr = (4S,) 2 + ipfXf- a©) 

When one more particle is extracted from the right side of Equation (j2J), we have 



S' 



where 



a nsm _ JN} {n-i} 

aNSM _ {N} h {N-l} ,{N} {N-l} [io\ 

A IJ,v,S ~ a SMn u S+l,M-^,f ^ u SMfj, a S-l,M-w> \ LO ) 

^^v,S-2 — °SMn°S~l,M-W \ lq: ) 



or A^Jf = otherwise. 
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Then the z-th state can be rewritten as 

e iM = J2x^)xu(2)^. (15) 

From the normality as before, we have 1 = X^i^^W Iv 9 ^ )• Thus, it is straight forward 
to define the 2-body spin correlative probability as 

= «K) = E dTdf,A»™A™;™. (16) 

SS'S" 

P l ^f is the probability that the spin of a particle is in [i while another in v when the 
two are observed, Obviously, P^f = Pl'jf . If more particles are extracted, higher order 
spin correlative probabilities could also be similarly defined. These probabilities do not 
have a counterpart in the MFT, therefore additional information might be provided by 
them. Incidentally, the technique for the measurement of the correlative probabilities 
is mature in particle physics and nuclear physics, but not in condensed matter physics. 
The development of related technique is desired. 

In general, one can define the correlation coefficient 7^ = P^/P^' P^' to 
measure quantitatively how large the correlation is. If 7^ deviates remarkably from 1, 
the correlation is strong. Whereas if j 1 '^ ~ 1, the correlation is weak, and the system 
can be well understood simply from the 1-body probabilities. 

Numerical examples will be given below. Tiu and \J%/ {muj) are used as units 
for energy and length, respectively, where m is the mass of atom. The spatial wave 
function <p(r) is obtained under the Thomas-Fermi approximation, and thereby we 
have the strength C2 = 0.154c2/(-/Vc ) 3 ^ 5 - Since a slight inaccuracy that may exist 
in 4>(r) would cause only a slight deviation in the magnitude of C2, the approximation 
is acceptable in the qualitative sense, u = 300 x 2ir (in sec -1 ) and N = 1000 are in 
general assumed (unless particularly specified). 87 Rb and 23 Na condensates will be used 
as examples for the C2 < and C2 > species, respectively. In the units adopted, we 
have c = 2.49 x lO -3 ^ and c 2 = -1.16 x lO -5 ^ for 87 Rb, and c = 6.77 x lO -4 ^ 
and C2 = 2.12 x lO^-^/u; for 23 Na. The diagonalization of -f^ mo d is straight forward 
when all the parameters are given. Then, the coefficients d 1 ^ 1 can be known. From 
Equations ([9]), ( !T0|) and ( Tl6|) . the 1-body and 2-body probabilities can be obtained. 

4. Low temperature limit 

The condensate will fall into its ground state Q\m with a specified M when T = 0. M 
is determined by how the species is prepared. The case with T ^ will be considered 
later. 

4-1. Condensates with Ci < 

To evaluate how strong the correlation in Oia/ would be, the correlation coefficients 
7o'o f of Rb against B are shown in Figure [H(a). The magnitudes of all the curves are 
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Figure 1. (Color online) The variation of the correlation coefficient Jq q 1 against B 
for the ground states of condensates with a specified M. N — 1000, w — 300 x 2tt 
are assumed, and the experimental data for 87 Rb and 23 Na are adopted (also for the 
following figures, except particularly specified). M is given at four values marked 
beside the curves. 



extremely close to one in (a). This is also true for r y 1 ^ 1 with (fi,v) ^ (0,0). It implies 
that the correlation in Q\m is very weak. 

When B = 0, Q\m will have S = N.^\ [15] Its expression does not depend on the 
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Figure 2. (Color online) P^' M of the 87 Rb condensate against B. (a) is for M = 
and fi is marked beside the curves, (b) is for [i — and M is marked beside the curves. 



details of interaction but only on the sign of C2, and can be uniquely written as 

M/ = $ N NM = (... ((x(l)x(2)) 2X (3)) 3 • • • x(N))n,m, (17) 
where (xaXb)x implies that the spins of xa and xb are coupled to spin A, and so on. 
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The special way of spin coupling in Equation ( |17p (i.e., the combined spin of an arbitrary 
group of j particles is j) assure that the state is normalized and symmetrized. Where 
all the spins are roughly aligned along a common direction, but the azimuthal angle of 
this direction is arbitrary. 

From Equation ( FlOl) the analytical forms of P^' M can be derived as 

, M (N±M)(N±M-1) 



P 



±1 



2N(2N- 1) 



nltM (N-M)(N + M) 

P ° ~ N(2N-1) • (19) 

Equation ( Fl9l) can be approximately rewritten as Pq' M ~ |(1 — 7^)- Thus, if M = 0, 
half of the particles would have fi = 0. If M ^ 0, the increase of M would lead to a 
decrease of Pq' M ■ In particular, Pq' M —> when M — > N as expected. Incidentally, the 
above analytical forms of P^ ,M are identical to those from the MFT when N — > oo.[6] 

When B varies, the one-body probabilities P^' M of Q\m against B are calculated 
and shown in Figure El The left ends of the curves coincide exactly with those from 
the above analytical expressions. Note that the fi ^ particles will gain additional 
energy from the quadratic Zeeman term. In order to reduce the energy the number of 
fi = particles will increase with B. This is clearly shown in Figure |2j For M = 
as shown in Figure Ett, the derivative of P 1,0 varies very swift from a positive value 
to zero at the vicinity of a turning point B = 0.28G where Pq 1 ' ~ 1, i,e, the number 
of \x = particles N is very close to its limit N. The swift variation at the turning 
point will become a sudden transition when N is large, and the point becomes a critical 
point P cr it. Passing through this point, the derivative varies abruptly, but the spin-state 
610 itself varies continuously. When B > B crit) the ground state remains unchanged 
as 0io = |0,iV,0) = (xo) N , where |JVi, JV , N -i) is a Fock-state with N x , N , and iV_i 
particles in \i = 1, 0, and -1, respectively. 

When M ^ 0, not all the particles can be changed to \i = because of the 
conservation of M. Accordingly, the increase of Pq' M is hindered and the curves in 
Figure [2b are smoother. When M is large, the curves are very flat implying that Qim is 
affected by B very weakly. When B — > oo, we found that 0i,m - > \M,N — M, 0). 
Accordingly, P 1,M — > (JV — M)/N, Pl' M -> M/N, and Plf ->■ 0. In this way, 
N = N — M is maximized so that the quadratic Zeeman energy is minimized. 

In conclusion of this section for small condensates of Rb, the ground state 6i m is 
continuously changed from t° N — M, 0) without a transition when B increases. 
However, the derivative undergoes a transition at B crit when M — 0. It is recalled 
that the MFT, which is correct when N is very large, has predicted the transitions 
between the ferromagnetic, broken-axisymmetry, and the polar phases (3J, [HI E]- The 
broken-axisymmetry phase is caused by a rapid quenched field, therefore it is not 
expected to appear in our case with B varying adiabatically.. 
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4-2. Condensates with c 2 > 

Contrary to the previous case, the spin correlation is very strong in the ground states 
0i,M of C2 > condensates as shown by ^q^ 1 in Figure [lb, where the left ends of 
all the curves are remarkably larger than 1. It implies that, when two particles are 
observed, the probability of a particle in fi = would be much larger if the other one 
has v — 0. Starting from B = there is a domain of B where 7o'o^ remains large and 
nearly unchanged (refer to the M = 100 curve in Figure [lb). Outside the domain, the 
correlation vanish rapidly. The domain is called a strong correlation domain (SCD), 
where the spin- structure at B = will keep unchanged against the increase of B. The 
SCD does not have a clear border. For each curve in Figure [Tb, let _B sc d be the value 
where 7 ' J = [( 7o ;o )max + l]/2. Then, B scd is considered as the outward border of the 
SCD. A larger M will lead to a larger B scd (say, for Figure [lb, if M = 10, 100, and 500, 
respectively, _B scd 0.0041G, 0.0318G, and 0.153G). Thus the spin structure at B = 
will better keep when M is large. 

When B = 0, the ground state will have S as small as possible. [5j [15] Thus, 
©i,m = $m m H N — M is even, or = "&m+xm ^ N — M is odd. Due to the fact 
that $sm i s unique, we have 

KM^n(Xi) M [(xxU iN - M)/2 h (20) 

where V is the symmetrizer, i.e., a summation over the N\ particle permutations. So 
this state is composed of a group of \i = 1 particles together with a group of singlet 
pairs. And 

K + i,m oc n«U[(xx)oF- M - 1)/2 }, (21) 
which is composed of a group of M + 1 polarized particles (but the direction of 
polarization deviates a little from the Z-axis) together with a group of singlet pairs. 
From Equation ffTOj) the one-body probabilities are 



p i,m = |M|(2 + l/AQ-2|M| 2 /iV -1 

" (2|M| +3)(2|M| - 1) ' 1 ' 

p x,M = |M|(2 - 1/AQ + |M| 2 (4 + 2/AQ ± M(4|M|(|M| + 1) - 3)/N -2 

±x " 2(2|M| +3)(2|M| - 1) ^ } 

For the case M = and iV being even, Q± : m = ^oo ^ 'P([(xx)o] N ^' 2 - This 
state is named the polar state where all the particles are in the singlet pairs. From 
Equations. ( 122]) and ([23"]) . the polar state has P^'° = 1/3 for all fi. In other words, 
the particles are equally populated among the three components. This is a common 
property of S = states. 

When B increases, P^' M against B are shown in Figure [3J Since the ground state 
will have more and more \x = particles, P ' increases sharply from 1/3 to ~ 1 as 
shown in Figure^. Correspondingly, the ground state B^o is sharply transformed from 

to \0,N, 0) = {xo) N ■ Comparing the curve of Pq'° in Figure [3^l with the one in 
Figure [2^ for Rb, the rise of the former is much faster than the latter. Thus the Na 
condensate is highly sensitive to the appearance of B if M = 0. A very weak field (a 
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few mG) is sufficient to break nearly all the pairs and turn every spin lying on the X-Y 
plane independently. The independence is supported by the curve of 7 ' in Figure [Tb, 
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where 7 g ~ 1 except B is close to zero. However, the high sensitivity will be lost when 
M^O. 
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Figure 4. (Color online) Correlative probabilities P^„ oa of the 23 Na condensate 
against B. (//, v) is marked beside the curves. 

In Figure |3b the curve of P 1,100 remains to be horizontal when B < B sc & = 0.0318G 
(refer also to the curve with M = 100 in Figure [lb). Thus the spin-structure inside 
the SCD might remain unchanged. To clarify, the 2-body probabilities are calculated 
and shown in Figure HJ where all the remain unchanged in the SCD. Thus the 

invariance of the structure against B is strongly supported. The invariance can also be 



seen by observing the overlap (e^°|9f M ) 2 = ($m,m\®i,m) 2 = (4/ f against B (N 
and M are assumed to be even, a superscript B is added to emphasize the dependence 
on B). With the parameters of Figure [3] and with M = 100, this overlap is > 0.99 
when B < 0.0267. Thus the invariance is directly confirmed. However, it decreases very 
fast when B is close to -B sc d, and is equal to 0.58 when B = -B SC( j. Afterward, it tends 
to zero rapidly when B is further larger. The existence of the SCD demonstrates that 
the mixture of a group of singlet pairs together with a group of M unpaired particles 
(each has fi = 1) is capable to keep its structure against B. However, the capability 
will be lost when B is close to or > -B sc d- Afterward the mixture will begin to change. 
Note that the change is characterized by the increasing of [i = particles, which come 
from the breaking of pairs. Therefore the change would be less probable if the original 
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number of pairs is small. Thus a larger M (implying a fewer original pairs) will lead 
to a better stability and therefore a larger P sc d- A larger c 2 will also lead to a better 
stability and therefore a larger P sc d (say, if c 2 is one time larger than the experimental 
value of (c 2 )Na, then the P sc d would be enlarged from 0.0318G to 0.0465G for the curve 
with M = 100 in Figure [3b ). 

Since M affects the stability, P^' M is in general sensitive to M. If B is weak, the 
sensitivity would be very high when M is small. E.g., it is shown in Figure |3b that Po' M 
will decrease dramatically at B = 5mG when M/N is simply changed from to 0.01. 

When B — > oo, all the pairs will be destroyed and ©i 5 a/ will tend to \M, N — M, 0). 
Thus, disregarding c 2 < or > 0, both species tend to the same structure. 

4-3. Fidelity susceptibility 

In order to understand the sensitivity of the ground states Of M against B quantitatively, 
the fidelity susceptibility [T6| [T71 [18] . 

r M (B) = iim^(i - |(ef+/|ef M )|), (24) 

is calculated and given in Figure [5] This figure demonstrates that each species has its 
own region highly sensitive to B. For Rb, the most sensitive region is surrounding the 
critical point P cr i t where the increase of N stops suddenly (refer to Figure Eh). When B 
is small, the sensitivity does not depend on M and is in general very weak. Whereas for 
Na, the most sensitive region is surrounding the border of SCD, P sc d, where Nq begin 
to increase suddenly (refer to Figure [3b). The fact that P sc d will become larger with M 
is clearly shown in Figure [5b. The sensitivity can be very high when B is small if M is 
also small. 

4-4- Effect of the trap and the particle number 

Since the linear Zeeman term does not affect the spin-structures, the ratio g/C 2 involved 
in H mo d is crucial. This quantity is proportional to B 2 /ui 6 ^ 5 . Therefore, a larger to would 
reduce the effect of B. Consequently, for oj' > to, all the curves plotted in Figure [1] 
to Figure H] will extend horizontally to the right by a common factor (tu'/tu) 3 ^ 5 . In 
particular, the SCD will become larger (say, B scd of the ground states O^ioo of Na 
would increase from 0.0318G to 0.0655G if u is from 300 x 2vr to 1000 x 2vr). 

On the other hand, the change of iV causes not simply a change of scale, because the 
numbers of degrees of freedom are thereby changed. For an example, if N = 1000, 100, 
and 10, P 1,0 of Rb would be = 0.99 when B = 0.277G, 0.178G, and 0.134G, respectively. 
Thus, the curve of P 1,0 in Fi gure [2k would rise up faster if iV decreases from 1000. On 
the contrary, if iV = 1000, 100, and 10, P 1,0 of Na would be = 0.90 when B = 0.003G, 
0.018G, and 0.070G, respectively. Thus, the curve of P 1,0 in Fi gure [3^ would rise up 
slower if iV decreases from 1000. Thus, the effect of iV on c 2 < and > species is 
different. 
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Figure 5. (Color online) ln{T m{B)) against M/N and B for the ground states of 
Rb (a) and Na (b) condensates. B is in Gauss and N — 1000. The maximum of 
ln{T m{B)) = 10.65 (a) or 11.20 (b). The region surrounding the maximum is marked 
with "X". The inmost contour of the X-region has the value 9.275 (a) or 9.078 (b). 
The difference of the values of a contour and its adjacent outer contour is 1.375 (a) or 
2.122 (b). 

The effect of N on Pq' M with M ^ is shown in Figure El The curve with N = 1000 
in Figure [6^ is identical to the curve with M = 100 in Figure [2b for C2 < 0. In Figure^ 
for Rb the curves with a smaller N will rise up faster against B just as the previous case 
with M — 0. For C2 > 0, the curve with N = 1000 in Figure [6b is identical to the curve 
with M = 100 in Figure [3b. Similar to the case with M = 0, a smaller iV will cause 
also a smoother change. In particular, the abrupt change appearing in the vicinity of 
5 sc d as shown in Figure [3b will disappear when iV is small. 

5. Finite temperature 

Since the level density in a condensate is usually dense, thermo-fluctuation is in general 
not negligible. What actually measured is the weighted probabilities Pjf = J2i WiP l ^ M 
and = Ei WiPjtf, where W { = ex P (-E t /k B T)/ e W {-E 3 /k B T) is the weight, E t 
the energy, T the temperature, and the summation in principle runs over all the states. 
However, when T << hco/ks = T , the contribution arises essentially from the ground 
band, and all the higher states can be neglected. [15J The members of the ground band 
have nearly the same spatial wave functions, and their spin-states together with Ei can 
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Figure 6. (Color online) Pq' m against B with M = N/10 and N is given at five values: 
N = 1000(solid line), 400(dash), lOO(dash-dot), 40(dash-dot-dot) and 10 (solid). 



be obtained via the diagonalization of H mo d- 

As an example, Pq ^ of the condensates with C2 < at T = To/10 is plotted in 
Figure Uh- In order to see the effect of interaction, is given at three values. When B 
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(b) c. <0 and B = 0.2G 
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Figure 7. (Color online) The weighted probabilities for the condensates with C2 < 0. 
T = To/10 is assumed and cq = (co)Rb (i.e., the experimental value of 87 Rb) and 
C2 = Qi(c2)ab adopted, (a) P against _B with a given at three values marked 
beside the curves. The dash curve has a — 1. (b) Pq '' and P q against a with 
S = 0.2G. 



is very small (< 0.0267) or sufficiently large (> 0.35G), the three curves associated with 
different values of C2 overlap nearly. However, in between, they are rather sensitive to c 
as shown in Figure [7^.. To see clearer, both Pq ^ and against C2 under B = 0.2G are 
plotted in Figure {7b- The figure demonstrates that the weighted 2-body probabilities are 
very sensitive to Ci- Furthermore, the patterns of the curves can be tuned by altering B. 
Therefore, the measurement of the weighted probabilities under various B can provide 
rich information on the parameters of interaction. 

One more example for the case C2 > is shown in Figure El Similar to the previous 
case, high sensitivity to the interaction is found. In Figure [Sk, the left ends of the curves 
are flat. This is caused by the existence of the SCD, where the spin- structure remains 
nearly unchanged. 

6. Final remarks 

We have studied the spin structures of small spin-1 condensates (N < 1000) under a 
magnetic field B. The theory is beyond the MFT, and the single mode approximation 
has been adopted. The fractional parentage coefficients have been used as a tool for the 




Figure 8. (Color online) Similar to Figure [7] but for the condensates with C2 > and 
M = 100. c = (c )Na and c-i = a(c2)Na are adopted. Refer to Figure [7] 



analytical derivation. The 1-body, 2-body, and weighted probabilities are defined and 
calculated to extract information from the spin eigenstates. The correlation coefficients 
7^f J and the fidelity susceptibility Tm(B) have also been calculated. The following 
results are mentioned (N — M is assumed to be even for simplicity). 

(i) When B = 0, the ground state Qi t M is either the ferromagnetic state fi^M ^ 
C2 < 0, or is the polar state M if c 2 > 0. For fi^Mi a ^ the N spins are coupled to total 
spin S = N, 7^ M is close to 1 and r M (0) is very small (Figured] and Figure [5]), therefore 
the ferromagnetic state does not have spin-correlation, and is inert to B when B is small. 
The state ~&mm ^ as M particles in xi together with (N — M)/2 singlet pairs, its •j^q 1 
is very large. Thus the polar state contains strong spin-correlation. Furthermore, its 
Tm{B = 0) is very large when M is small. It implies a high sensitivity against the 
appearance of B. 

(ii) When B increases, the number of /i = particles Nq in Qx,M wm i n general 
increase so as to reduce the quadratic Zeeman energy. For c 2 < 0, Oi,m will be changed 
from t° \M, N — M, 0) when B is from to oo. The change goes on continuously, 
no transition in spin-structure occur. In accord with the change of 0i,m, No is changed 
from (Jy ~ffy M) ~ (1+M/N) ^ N _ M ) ( refer to Equation flH])) to N — M. Therefore, 
(No) b^qo-(No) b=o ~ (N—M) 2 /2N = N 0) diff, which is the maximal number of particles 
allowed to be changed from being \i =fi to // = 0. When M is close to N, N 0t diff is 
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very small implying that the room left for changing is very small, therefore Oi 5 m is 
inert to B (refer to Figured where the curve with M = 980 is very flat). When M is 
smaller, Gi,m has much room for changing and therefore would be more sensitive to B. 
In particular, when M = 0, there is a critical point B crit . Once B > B crit (No)b — N 
and the ground state varies with B no more. There is an abrupt change in the derivative 
at -Bcrit- Thus, although the spin-structures vary continuously with B, the related 
derivatives might not. B ciit is equal to 0.28G in Figure [2j It will be smaller when N 
decreases, and larger when u increases. 

(iii) For c 2 > , the increase of B from to oo causes a change of 0i,m from T?j^ M 
(or M ) to \M, N — M, 0). 0i,m of both cases c 2 < and > tend to the same state 
because it is the most advantageous state for reducing the quadratic Zeeman energy. The 
change of 0i,m goes on also continuously without transitions. In this process N is in 
general increasing via a mechanism, i.e., a breaking of pairs as (xx)o XoXo (whereas 
the process (xx)o ~>* XiX-i is suppressed under B because it causes an increase of 
quadratic Zeeman energy). It is recalled that the particle correlation in 0i,m is very 
weak when c 2 < 0, but strong when c 2 > due to the formation of pairs. When all the 
particles are paired (i.e., M = 0), the structure is extremely sensitive to B. A very weak 
B (a few mG) is sufficient to break all the pairs. For a comparison, Pq 1 ' in Figure [2^ 
for Rb will be equal to 0.9 when B = 0.25G, but only = 0.003G in Figure [3^ for Na. 
However, when unpaired particles emerge (i.e., M^0), the pairs will have an additional 
ability to keep themselves. The mechanism underlying this phenomenon deserves to be 
studied further. This leads to the appearance of SCD ranging from B = to -B sc d- A 
larger M will lead to a larger -B scc i, while 5 scc i = when M = 0. The spin-structure 
will remain unchanged when B < B sc a, but Nq will begin to increase when B > B SC( i. 
The derivative of the probabilities varies very swiftly in the vicinity of -B sc d- The swift 
variation will become a sudden jump when iV is large (refer to Figure |6b). Thus _B sc d 
is also a critical point when iV is large. As a numerical example, when iV = 1000 and 
M = 100, -B SCC ! = 0.03G in Figure |3b. Similar to B CTit , B sc ^ will become smaller when TV 
decreases (with M/N remaining unchanged), and will become larger when u increases. 

(iv) The sensitivity of the ground states against B is quantitatively shown in 
Figure 

(v) When B is appropriately chosen and T is sufficiently low, the measurable 
weighted probabilities may provide rich information on the parameters of interaction. 
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